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AND RIEMANN-HILBERT PROBLEM 



Alexei Borodin 

Abstract. We show how the Riemann-Hilbert problem can be used to compute 
correlation kernels for determinantal point processes arising in different models of 
asymptotic combinatorics and representation theory. The Whittaker kernel and the 
discrete Bessel kernel are computed as examples. 



Introduction 

A (discrete or continuous) random point process is called determinantal if its 
correlation functions have the form 



= 15 



where K{x,y) is a function in two variables called the correlation kernel. A major 
source of such point processes is Random Matrix Theory. All the "unitary" or 
"/3 = 2" ensembles of random matrices lead to determinantal point processes which 
describe the eigenvalues of these matrices. 

Determinantal point processes also arise naturally in problems of asymptotic 
combinatorics and asymptotic representation theory, see [B01]-[B04], [BOO], [J], 
[012]. Usually, it is not very hard to see that the process that we are interested 
in is determinantal. A harder problem is to compute the correlation kernel of this 
process explicitly. The goal of this paper is to give an informal introduction to 
a new method of obtaining explicit formulas for correlation kernels. It should be 
emphasized that in representation theoretic models which we consider the kernels 
cannot be expressed through orthogonal polynomials, as it often happens in random 
matrix models. That is why we had to invent something different. 

The heart of the method is the Riemann-Hilbert problem (RHP, for short). This 
is a classical problem which consists of factorizing a matrix-valued function on a 
contour in the complex plane into a product of a function which is holomorphic 
inside the contour and a function which is holomorphic outside the contour. It 
turns out that the problem of computing the correlation kernels can be reduced to 
solving a RHP of a rather special form. The input of the RHP (the function to be 
factorized) is always rather simple and can be read off the representation theoretic 
quantities such as dimensions of irreducible representations of the corresponding 
groups. We also employ a discrete analog of RHP described in [B2]. 

The special form of our concrete RHPs allows us to reduce them to certain linear 
ordinary differential equations (this is the key step), which have classical special 
functions as their solutions. This immediately leads to explicit formulas for the 
needed correlation kernels. 
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The approach also happens to be very effective for the derivation of (nonhnear 
ordinary differential) Painleve equations describing the "gap probabilities" in both 
random matrix and representation theoretic models, see [BD], [B3]. However, this 
issue will not be addressed in this paper. 

The paper is organized as follows. In §1 we explain what a determinantal point 
process is and give a couple of examples. In §2 we argue that in many models 
correlation kernels give rise to what is called "integrable integral operators" . In §3 
we relate integrable operators to RHP. In §4 we derive the Whittaker kernel arising 
in a problem of harmonic analysis on the infinite symmetric group. In §5 we derive 
the discrete Bessel kernel associated with the poissonized Plancherel measures on 
symmetric groups. 

This paper is an expanded text of lectures the author gave at the NATO Ad- 
vanced Study Institute "Asymptotic combinatorics with applications to mathemat- 
ical physics" in July 2001 in St. Petersburg. It is a great pleasure to thank the 
organizers for the invitation and for the warm hospitality. The author would also 
like to thank Grigori Olshanski and Percy Deift for helpful discussions. 

This research was partially conducted during the period the author served as a 
Clay Mathematics Institute Long-Term Prize Fellow. This work was also partially 
supported by the NSF grant DMS-9729992. 

1. Determinantal point processes 

Definition 1.1. Let X be a discrete space. A probability measure on 2^ is called 
a determinantal point process if there exists a function K : X x X ^ C such that 

Prob{^ e 2^ 1^ D {xi, . ..,Xn}} = det[K{xi,Xj)]lj^-^ 

for any finite subset {xi, . . . .Xn} of X. The function K is called the correlation 
kernel. The functions 

Pn : {n-point subsets of X} — > [0, 1] 
pn : {xi,.. .,Xn} ^ Prob{^ 1^ D {xi,...,Xn}} 

are called the correlation functions. 

Example 1.2. Consider a kernel L : X x X ^ C such that 

• det[L(xi, Xj)]'lj^i > for all fc-point subsets {yi, . . . , i/k} of X. 

• L defines a trace class operator in for example, J2x yex l-^i^^v)] < 
or L is finite rank. In particular, this condition is empty if \X\ < oo. 

Set 

Prob {{yi, . . . , yk}} = ^^^^^ ^ • det[L(yi, yj)]lj=i . 

This defines a probability measure on 2^ concentrated on finite subsets. Moreover, 
this defines a determinantal point process. The correlation kernel K{x,y) is equal 
to the matrix of the operator K = L{1 + L)-^ acting on e'^{X). See [DVJ], [BOO, 
Appendix] for details. 
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Definition 1.3. Let X be a finite or infinite interval inside R (e.g., R itself). A 
probability measure on locally finite subsets of X is called a determinantal point 
process if there exists a function K : X x X ^ C such that 



^. ProbjA e 2^^ I A intersects [xi, Xi + Axi] for alH = 1, . . . , n} 

Aa;i,...,Aa;„^0 Axi • • • Ax^ 

= Ael[K{xi,Xj)]l^^^ 

for any finite subset {xi, . . . ,a;„} of X. The function K is called the correlation 
kernel and the left-hand side of the equality above is called the nth correlation 
function. 

Example 1.4. Let w{x) be a positive function on X such that all the moments 
J^x'^w{x)dx are finite. Pick a number N and define a probability measure on 
iV-point subsets of X by the formula 

PN{dxi, . . . ,dxN) = cn Y\. {xi-XjY W w{xk)dxk. 

l<i<j<N l<k<N 

Here cjv > is a normalizing constant. This is a determinantal point process. 
The correlation kernel is equal to the A'"th Christoffel-Darboux kernel KN{x,y) 
associated with w{x), multiplied by ^/w(x)w(y). That is, let 

po = 1, pi{x), p2{x),... 

be monic (= leading coefficient 1) orthogonal polynomials on X with the weight 
function w{x): 

Pmix) = x^ + lower order terms , 

/ Pm{x)Pn{x)w{x)dx ^ hmSmn, m, n = 0, 1, 2, . . . . 

Jx 

Then the correlation kernel is equal to 



K.{x,y) = j:'-^^^V^IK^ 
fe=o 

^ 1 PN{x)pN-l{y) -pN-i{x)pN{y) 
hN-1 x-y 



\/w{x)w{y). 



The construction of this example also makes sense in the discrete setting. See [Dy] , 
[Me], [NW], [J] for details. 

Remark 1.5. The correlation kernel of a determinantal point process is not defined 
uniquely! In particular, transformations of the form K[x,y) — > j^K(x,y) do not 
change the correlation functions. 
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2. Correlation kernels as integrable operators 



Observe that the kernel Kn[x, y) of Example 1.4 has the form 

(j){x)il;{y) - ip{x)(f){y) 



KN{x,y) = 



x-y 



for appropriate cf) and -0. Most kernels appearing in "/9 = 2 ensembles" of Ran- 
dom Matrix Theory have this form, because they are either kernels of Christoffel- 
Darboux type as in Example 1.4 above, or scaling limits of such kernels. However, 
it is an experimental fact that integral operators with such kernels appear in many 
different areas of mathematics, see [De]. 



Definition 2.1. An integral operator with kernel of the form 

h{x)gi{v) + f2{x)g2{y) 
x-y 



(2.1) 



is called integrable. Here we assume that fi{x)gi{x) + f2{x)g2{x) = so that there 
is no singularity on the diagonal. Diagonal values of the kernel are then defined by 
continuity. 

The class of integrable operators was singled out in the work of Its, Izergin, 
Korepin, and Slavnov on quantum inverse scattering method in 1990 [IIKS]. 

We will also call an operator acting in the £^-space on a discrete space integrable 
if its matrix has the form (2.1). It is not obvious how to define the diagonal entries 
of a discrete integrable operator in general. However, in all concrete situations we 
are aware of, this question has a natural answer. 

Example 2.2 (poissonized Plancherel measure, cf. [BOO]). Consider the 
probability measure on the set of all Young diagrams given by the formula 



Prob{A} = e-'^^l-'l (^^^ J . (2.2) 

Here ^ > is a parameter, dim A is the number of standard Young tableaux of 
shape A or the dimension of the irreducible representation of the symmetric group 
S\\\ corresponding to A. Denote by (pi, . . . ,^^191, . . . ,qd) the Frobenius coordinates 
of A (see [Ma, §1] for the definition of Frobenius coordinates). Here d is the number 
of diagonal boxes in A. Set Z' = Z + | = {=t|; =t|; ■ ■ ■ }• 

Let us associate to any Young diagram A = (p | g) a point configuration Fr(A) C 
Z' as follows: 

= {Pi + • • • + |, -Q-i - |, • • • , -9d - 1} • 

It turns out that together with (2.2) this defines a determinantal point process 
on Z'. Indeed, the well-known hook formula for dim A easily implies 



Prob{A} = e"^ det 



2 

e— 



ip^-mQJ-mp^+qJ) 



= e-'det[L{y,,yj)]^A, 



where {j/i, . . . , j/2d} = f^l^^)) ^'Hd L{x,y) is a Z' x Z' matrix defined by 



0, 



L{x,y) 



„M±M 
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if xy > 0, 
, if xy < 0. 



In the block form corresponding to the sphtting Z' = Z^ U Z'_ it looks as follows 



L{x,y) 



9' 



-x+y 



-X 



my-l)lx-y 



The kernel L{x,y) belongs to the class of integrable kernels. Indeed, if we set 



fi{x) ^g2{y) 



0, 



a;<o, I {-x-\y: 



a; > 0, 
a; < 0, 



then it is immediately verified that L{x,y) = {fi{x)gi{y) + f2{x)g2{y))/ {x — y). 
Comparing the formulas with Example 1.2, we also conclude that = det(l + L)} 

What we see in this example is that L is an integrable kernel. We also know, 
see Example 1.2, that the correlation kernel K is given hy K = L(l + L)~^ . Is this 
kernel also integrable? The answer is positive; the general claim in the continuous 
case was proved in [IIKS] , the discrete case was worked out in [B2] . 

Furthermore, it turns out that in many situations there is an algorithm of 
computing the correlation kernel if L is an integrable kernel which is "simple 
enough" . The algorithm is based on a classical problem of complex analysis called 
the Riemann-Hilbert problem (RHP, for short). 

Let us point out that our algorithm is not applicable to deriving correlation 
kernels in the "/3 = 2" model of Random Matrix Theory. Indeed, the Christoffel- 
Darboux kernels have norm 1, since they arc just projection operators. Thus, it is 
impossible to define the kernel L = K{1 — K)~^ , because (1 — K) is not invertible. 
In this sense, RMT deals with "degenerate" determinantal point processes. 

On the other hand, the orthogonal polynomial method of computing the corre- 
lation kernels, which has been so successful in RMT, cannot be applied directly to 
the representation theoretic models like Example 2.2 above (see, however, [J]). The 
algorithm explained below may be viewed as a substitute for this method. 



3. Riemann-Hilbert problem 

Let E be an oriented contour in C. We agree that (-l-)-side is on the left of the 
contour, and (— )-side is on the right of the contour. Let v be a 2 x 2-matrix valued 
function on S. 



^Since '^^ix j/ez' l-^(*>l')l < oOj the operator L is trace class, and det(l + L) is well-defined. 
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Definition 3.1. We say that a matrix function m : C \ E — > Mat(2, C) solves the 
RHP (E, v) if 

(1) m is analytic in C \ E; 

(2) 171+ = rri-V on E, where 'm±{x) = Um "^(0- 

»a; from (±)-side 

We say that m solves the normalized RHP (E, v) if, in addition, we have 

1 0' . 
as C 



(3) m(C) 



1 



oo. 



Next we explain what is a discrete Riemann-Hilbert problem (DRHP, for short). 
Let X be a locally finite subset of C, and let w be a 2 x 2-matrix valued function 
on X. 

Definition 3.2. We say that a matrix function m : C \ X — > Mat(2, C) solves the 
DRHP (X, w) if 

(1) m is analytic in C \ X; 

(2) m has simple poles at the points of X, and 

Kes^^x'iTT'iO — lini(m(C)ty(a;)) for any x E X. 

We say that m solves the normalized DRHP (X, w) if 

1 0^ . 

as C — OO- 



(3) m(C) 



1 



If the set X is infinite, the last relation should hold when the distance from ( to 
X is bounded away from zero. 

Our next step is to explain how to reduce, for an integrable operator L, the 
computation of the operator K = L(l + L)~^ to a (discrete or continuous) RHP. 

Continuous picture [HKS]. Let L be an integrable operator on -L^(E, |(iC|), E C 
C, with the kernel {x, j/ G E) 



L{x,y) = 



fi{x)9i{y) + f2{x)g2{y) 



fiix)giix) + f2{x)g2{x) = 0. 



x-y 

Assume that (1 + L) is invertible. 

Theorem 3.3. There exists a unique solution of the normalized RHP (E, v) with 



V = I + 2ni 
For a; e E set 



7i" 






[gi 92] = 



1 + 2TTifigi 27iifig2 
27ri/2fi'i 1 + 2mf2g2 



Fi{x) 
F2{x) 

Gi{x) 
G2{x) 



- lim m{() 
lim m~*(C) 

C^x 



h{x) 

9i{x) 
92{x) 



Then the kernel of the operator K = L(l + L) ^ has the form {x, y E T,) 



K{x,y) = 



F,{x)G,{y) + F2{x)G2{y) 



x-y 



and Fi{x)Gi{x) + F2{x)G2{x) = 0. 



6 



Example 3.4. Let E be a simple closed curve in C oriented clockwise (so that the 
(+)-side is outside E), and let L be an integrable operator such that the functions 
/17 12^91^92 can be extended to analytic functions inside E. Then the solution of 
the normalized RHP (E, v) has the form 



m 



1 
1 

I -27ri 



fi 

/2 



outside E, 
[91 92] inside E. 



Then we immediately obtain Fi = /j, Gi = gi, i = 1, 2; and K = L(l + L)~^ = L. 
On the other hand, this is obvious because Jj,L{x,y)L{y,z)dy = by Cauchy's 
theorem which means that = 0. 

Discrete picture [B2]. Let L be an integrable operator on £^{X), X C C, with 
the kernel 

fi{x)9i{y) + f2{x) g2{y) , 

1 X ^ y, 



L(x,y) = 



x-y 



x = y, 

with fi{x)gi{x) + f2{x)g2{x) = 0. Assume that (1 + L) is invertible. 
Theorem 3.5. There exists a unique solution of the normalized DRHP {X, w) with 



w = — 



fi 

/2 



[91 92. 



-fi9i -fi92 

-f29l -f292 



For X & X set 



Fi{x) 
F2{x) 

Gi{x)' 
G2{x) 



= lim m{Q 
= lim m~*{Q 



fi{x 

f2{x 

9i{x) 
9i{x) 



Then the kernel of the operator K = L(l + L) ^ has the form {x, 2/ G E) 

F,{x)G,{y) + F2{x)G2{y) 



K{x,y) 



■< II 

[Gi{x) G'2(a;)]lim (m'(C) 



fi{x) 

Vf2{x) 



xi-y, 
x = y. 



We also have Fi{x)G\{x) + F2{x)G2{x) = on X . 

Theorem 3.5 can be extended to the case when L(a;, x) ^ 0, see [B2, Remark 
4.2]. 

Example 3.6. Let X = {a, b} be a two-point subset of C, and 



Then L is integrable with 




/i 



iy{b — a) 



/2 = 





1^ 



/i(a — b) 
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91 = 



92 = 



The notation means that, say, /i(a) = 0, fi{b) = u{b — a). Then 



w{a) 



(1 

l^{b-a) 



w{b) = 



u{a - b) 




Then the matrix m{Q has the form 

m{0 = I+- ^ ^ 

C — a 1 — //z/ 

One can check that detm = 1, and 

— a 1 — //i^ 



m-*(C) = / + 



1/ 
-1 



1 
J/ 



+ 



1 z/(6-a) 



C - 6 1 - //z^ 



+ 



1 iy{b-a) 



C-b 1- fxu 



-1 






1 




F2 = 



hm ( m'(^) 

C-»-o 



/i(a) 
/2(a) 




l—Hu 
1— /[ii/ 



hm ( m' iC) 
C^6 * 



/i(^>) 



l — fj,l/ 

1 

1— /it" 

,2 



l—flL' 
l—jJLV 



By Theorem 3.5, this implies that 

L 



K = 



1 



1 + L 1 - /xi/ 



V 



-jJLV 



which is immediately verified directly. Note that the condition 1 — jiv ^ is 
equivalent to the invertibility of (1 + L). 

In what follows we will demonstrate how to use Theorems 3.3 and 3.5 to compute 
correlation kernels of determinantal point processes arising in concrete representa- 
tion theoretic models. 

4. Harmonic ANALYSIS ON S'(oo): Whittaker kernel 

As is explained in [B02], see also [012], the problem of decomposing generalized 
regular representations of the infinite symmetric group 5(00) on irreducible ones 
reduces to computing correlation kernels of certain determinantal point processes. 

Specifically, consider a determinantal point process on Z' = Z + | constructed 
using Example 1.2 with the L- kernel given by 



L{x,y) 



0, 

\z{Z+l),_^i-Z+l)_y_i\C^ 

(a;-i)!(-y-|)!(a;-y) 

\zi-Z+l)_,_i{z + l)y_ijC 



xy > 0, 

a; > 0, y < 0, 

a; < 0, y > 0. 



Here z G C \ Z and ^ G (0, 1) are parameters. The symbol {a)k stands for a{a + 
1) • • • (a + /c - 1) = r(a + k)/T{a). 



Note that as |2;| — > oo, ^ — > 0, and \z\'^^ — > 6, this kernel converges to the L-kernel 
of Example 2.2. 

The problem consists in computing K = L{1 + L)~^ and taking the scaling limit 
]C{x, y) = lim(l - ■ K ([(1 - O-'x] + 1, [(1 - O'^y] + \) . 

This problem has been solved in [B02]. However, wc did not provide a derivation 
of the formula for the kernel K there, we just verified the equality K = L(l +L)~^. 

The goal of this section is to provide a derivation of the kernel /C(a;, y) bypassing 
the computation of K{x,y). 

Observe that there exists a limit 

Cix, y) = lim(l -i)-'-L ([(1 - i)-^x] + i, [(1 - i)-^y] + |) . 

Indeed, for a = [(1 - i)-^x], h=[{l- ^y^y], 
{z+l)a _ T{z + l + a) 



r(a + i) r(^ + i)r(a + i) r(z + i)' 
{-z^i)b _ r(-z + i + |6|) \h\-' 



r(6+i) r(-^ + i)r(|6| + i) r(-^ + i)' 

^§ ~ (1 - (1 - ^))2lT^ ~ e^, ^~^~e~2. 



x > 0, 

y<o. 



and we get 



0, 



sin7r2;| {x/\y\)'^^e 2 



-3: + y 



C{x,y) = < 



TT 2; y 

siuTT^I {y/\x\)^^e^^ 



TT 



x-y 



xy > 0, 
x>0, y <0, 

X <0,y> 0. 



It is natural to assume that /C = C{1 + C)~^. It turns out that this relation 
holds whenever £ defines a bounded operator in L^(M\ {0}), which happens when 
\^z\ < |, see [Oil]. Our goal is to derive /C using the relation C = jC{1 + jC)~^. 

It is easily seen that C is an integrable operator; we can take 



\z\ 2 



hix) =g2{x) 



f2{x) =gi{x) 



|r(^+i)l 



0, 



0, 



x^^e 2 ^ a; > 0, 



X <0, 
x> 0, 



\z\ 2 



\x\ 



"^^6 2 , X <0. 



\n-z+i] 

Note that C{y,x) = —C{x, y) which means that {1 + C) is invertible (provided that 
C is bounded). 

The RHP of Theorem 3.3 then has the jump matrix 



v{x) 



1 2i I sin TTZ I x^^^e 
1 

1 

2i\sm7rz\\x\~'^^^e^ 1 



x>0 
X <0. 



The key property of this RHP which will allow us to solve it, is that it can be 
reduced to a problem with a piecewise constant jump matrix. 
Let m be the solution of the normalized RHP (M \ {0}, v). Set 



*(C) = m(C) 







C~ 62 . 
Then the jump relation 171+ — m-V takes the form 



— Sffz 2. 





x'^'e 2 



f (x) 



S)f z — 2- 

x e 2 



X ■"^6 2 







and a direct computation shows that the jump matrix for takes the form 

a; > 



1 2i I sin TT^; 
1 

g27ri K2 Q 

2i|sin7r2;| e"^'^*^^ 



, a; < 0. 



Let us first find a solution of this RHP without imposing any asymptotic conditions 
at infinity. We will denote it by '^^ . Set 



*°(C), > 0, 

T 2z|sin7r2| 
1 



»C<o. 



Then has no jump across R-|-, and the jump matrix (\E'^) ^^P^]. on M_ has the 
form 



1 —2i\sim^z\ 
1 



2i\sm'Kz\ e"^"^*^^ 



+ 4|sin7r2;| —2i\snn:z\e' 



2i I sin 7TZ 



The determinant of this matrix is equal to 1, and the trace is eqiial to 2 cos(27r3?z) + 
4|sin7r2;p = e'^^^ + e~'^^^. Thus, if ^ ^ M, there exists a nondegenerate U such 
that 



+ 4|sin7r^| —2i\sm7rz\e 



-2-Ki Uz 



2i I sin T^z 



,-2Tvi ^z 



e-'^^^ 



U. 



This means that the matrix '^^U ^ has jump 



gTT^Z 

e"'^^^ 



across R_. Note 



that the matrix 



C'"^^ 



satisfies the same jump relation. Hence, 




is a solution of our RHP for ^. It foUows that ^'(^'°)~-'^ has no jump across M and 
this imphes, modulo some technicalities, that \E' = H"^^ where H is entire. 

Now we describe the crucial step. Since the jump matrix for \1/ is piecewise 
constant, = satisfies the same jump condition as and hence is 
meromorphic in C with a possible pole at C = 0. On the other hand we have 



1 
-1 



H ^ — entire function + 



A 



(4.1) 



where A has eigenvalues ±z Qz. 

Let us recall now that m solves the normalized RHP, which means that m(C) ~ / 
as \z\ — > oo. An additional argument shows that 



KC)=/ + mWC-^+0(|C|-^) 
with a constant matrix m^^\ Thus, 

*(c) = (/+mWr' + o (Id 

and 



1-2 



ICI 







oo, 



JRZg 2 



1 



*'(C)*-'(C) = -7^ 



"1 


" 


1 





-1 





^z 

(1) 



-m 



m. 



21 



(1) 

12 



+ o(|C|-'). 



Comparing this relation with (4.1) we conclude that 



1 
-1 



c 



Uz 
(1) 



m 



21 



(1) 
12 

-^Z 



*(C) 



*'(c) = (-^ ; + 

with mi'2m2i = (3?^)^ + (^C)^ = kP- This 1st order linear matrix differential 
equation leads to 2nd order linear differential equations on the matrix elements on 
, for example 



1 
2 




*ii- 



Using these differential equations and the asymptotics of \E' at infinity, it is easy 
to express in terms of the confluent hypergeometric function or the Whittaker 
function, see [Er, 6.9] for definitions. In terms of the Whittaker function W^^^^ the 
final formula for ^ has the form 



*(C) 



C--w^.+iMO \A (-C)-^W-sR.-i,.c,.(-C) 
-1^1 C--w^,-iMO (-C)-^^-iR.+i ^sj.(-C) 



It is not hard to show that det = 1, and 



-\z\{-o-'^w_ 



X-0 \z\C-^w^._r^,^AO 
(-0 C-^^^«.+i»3.(C) 
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Then Theorem 3.3 imphes 

1 

|r(i+i)| ^ii(^)^ 3;>0, 

1 

|r(-l+i)| ^i2(a;), x<0, 

1 

• |r(^+i)| ^21 (a^), a;>0, 



Fi{x) = 



|r(z+i)| ^2i(a^), a;>0, 

1 

|r(-l+i)| ^22(a:^), a:^<0, 

1 



Gi{x) 



and 



F2(X) - 



G2{X) 



|r(_^+-T)T*22(a;), a;<0, 

Fi(a;)Gi(|/) + F2(a:)G2(2/) 



IC{x,y) 



x-y 



x,yeR\{0}. 



This kernel is called the Whittaker kernel, see [BOl], [Bl]. 

5. POISSONIZED PLANCHEREL MEASURE: DISCRETE BESSEL KERNEL 

We now return to the situation described in Example 2.2. Our goal is to compute 
the correlation kernel K = L{1 + L)~^ . The exposition below follows [B2, §7]. 

According to Theorem 3.5, we have to find the unique solution of the normalized 
DRHP {Z\w) with 



w{x) — < 







g-x 



X e z'^, 
X e Z'_. 



Note that the kernel L is skew-symmetric, which means that (1 + L) is invertible. 
If we denote by m the solution of this DRHP then 



m(C) = / + 



a (3 
7 S 



C' + oi\c'\), ICI-00, 



with constant a, P, 7, S. The symmetry of the problem with respect to 



TOll(C) TOl2(C) 
TO2l(C) "7,22(0 



W22(-C) -"^21 (-C) 

-mi2(-C) "^ii(-C) 



implies that j — (3 and 6 — —a. 
Denote 77 = V9 and set 



niO = m(C) 





ry-^ 



Then n{Q solves a DRHP with the jump matrix 

1 




< 











12 



X e Z'+, 
, xeZ'_. 



12 



Note that this jump matrix does not depend on 77. This means that has the same 
jump matrix, and hence the matrix ^ is entire. Computing the asymptotics 
as C — > 00, we obtain 



dniO 
drj 



C -2/3 
2/3 -C 



+ 0(N-i). 



By Liouville's theorem, the remainder term must vanish, and thus 

dniO 



dr] 



2/3 -C 



n(C). 



(5.1) 



This yields 2nd order hnear differential equations on the matrix elements of n which 
involve, however, an unknown function (3 = f3{'q). 

In order to determine (5 we need to make one more step. Set 



P{0 = n{0 



r(C+i) 






1 



It is immediately verified that the fact that n solves the corresponding DRHP is 
equivalent to p being entire and satisfying the condition 



p{x) = (-1)^ ^p{x) 



1 

1 



X e Z'. 



(5.2) 



The key property of this relation is that it depends on x in an insubstantial way. 
This allows us to do the following trick which should be viewed as a substitute of 
the differentiation with respect to x. Set 



m = 



Pll(C + l) -Pi2(C + i) 

-P2l(C-l) P22(C-1) 



Then p satisfies the same condition (5.2) as p does, and hence 



r(C + i) 





n-c+h) 



satisfies the same DRHP as n does. Thus, nn ^ is entire. Working out the asymp- 
totics as C — > 00, we obtain 



n(C)n-nC) = 







4 



+oi\zr). 



Liouville's theorem implies that the remainder vanishes. Hence, 



Pll(C) =Pll(C + l) = ^P2l(C), P2l(C) = -P2i(C- 1) = -^Pii(C)- 
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This implies that {13/rf)'^ = 1- Both cases jS = ±77 lead, via (5.1), to linear 2nd 
order differential equations on the matrix elements of n or matrix elements of p. 
For example, (3 = —rj yields 



drf 



C(C-i) 



+ 4Un(C)=0, 



dvj 



General solutions of these equations can be written in terms of Bessel functions, 
and matching the asymptotics at infinity we obtain for fi — —rj 



and for /3 = ry 



p(C) = 



p(C) = 0? 



Jc-i(2r7) J-c+i(2r7) 



(5.3) 



Jc-i(2^) -J_c+i(2r7) 
^-C-i(2^) 



Here Ju{u) is the Bessel function, see [Er, 7.2] for the definition. 

Using the well-known relation = (— 1)'^J„ we immediately see that p{C,) 
given by (5.3) satisfies (5.2), while p(C) does not. In fact p{C,) satisfies (5.2) with 
(— 1)^~2 replaced with (— 1)^"'"2. This means that 



w(0 =p(0 







r?cr(-c + i: 

solves the initial DRHP, and by Theorem 3.5 we obtain 



Pii(x), X > 0, 

Pl2(x), X < 0, 

-P2l{x), X > 0, 

^22(2;), X < 0, 



F2{x) 
G2{x) 



P2l{x), X > 0, 
P22{x), X < 0, 

pii(x), X > 0, 

-Pl2(x), X < 0, 



and 



X - 2/ 



The diagonal values -ftr(x,x) are determined by the L'Hospital rule: 

K{x, x) = F[{x)Gi{x) + F^(x)G2(x). 

This is the discrete Bessel kernel obtained in [BOO]. The restriction of K{x,y) 
to X Z'_j_ was independently derived in [J]. 

It is worth noting that the matrix p also has an important meaning. In fact, 
if we define a kernel K using the formulas above with p replaced by p then K — 
L(L-l)-i. 
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